Formation of a Metallic Contact: Jump to Contact Revisited 
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The transition from tunneling to metallic contact between two surfaces does not always involve 
a jump, but can be smooth. We have observed that the configuration and material composition 
of the electrodes before contact largely determines the presence or absence of a jump. Moreover, 
when jumps are found preferential values of conductance have been identified. Through combination 
of experiments, molecular dynamics, and first-principles transport calculations these conductance 
values are identified with atomic contacts of either monomers, dimers or double-bond contacts. 
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Matter at the atomic scale has attracted much atten- 
tion during the last two decades. Not only because of 
the new properties arising when size is decreased, but 
also because the mechanisms found there can help us un- 
derstand behavior at the macroscopic scale. This is the 
case of atomic contacts where many new and unexpected 
phenomena have been found [l| and, at the same time, 
have provided clues about many technological problems 
in the macroscopic scale. A clear example would be the 
study of adhesion or friction which are consequences of 
the formation of a contact between two bodies, a pro- 
cess which always involves the formation of at least one 
atomic contact [2]. 

The formation of a contact between two metallic elec- 
trodes can be studied by measuring the electronic trans- 
port through them. It has been observed that right be- 
fore the electrodes are brought into contact, electrons can 
tunnel from one electrode to the other. This tunneling 
current, as shown by the curve on the left in Fig. HJ will 
increase exponentially while the electrodes come together 
and at some point a Jump to Contact (JC) occurs: the 
electrode separation becomes unstable and a sudden in- 
crease of the current takes place because of the formation 
of a single- atom contact [l|,|3| with a conductance of about 
a quantum unit(Go = 2e 2 /h). When the electrodes are 
then separated, a hysteretic behavior is observed in which 
the contact breaks in a jump to form a vacuum gap while 
relieving the elastic tension accumulated [4] . 
The process described above has been seen as the gen- 
eral rule in the formation of a metallic contact with a 
few exceptions to be understood [1]. It has been reported 
that in some cases the jump does not occur and, instead, 
there is a continuous increase of current from tunneling to 
contact, after which it saturates forming a plateau show- 
ing no hysteresis when the electrodes are separated. Ex- 
amples of this behavior have accumulated and presently 
have been reported for the cases of a W(lll)-Au(lll) 
contact [5], W nanocontactsQ, and for the case when 
a STM tip approaches Ag and Cu adatoms on (111) 



surfaces @- 

Here, the formation of a metallic contact for various met- 
als (Au, Pt, Ag, Ir, W, and Ni) is investigated. A sta- 
tistical analysis using density plots has been developed 
that, combined with atomistic simulations and ab-initio 
calculations of the electronic transport, provide new in- 
sight into the problem. The absence of a JC is not an 
exception but it occurs for many metals, in particular for 
Ni, W and Ir. 

In order to control the formation of the atomic contacts 
we have employed a home-made high stability Scanning 
Tunneling Microscope (STM) and Mechanically Con- 
trolled Break Junctions [1] (MCBJ). In STM experiments, 
the tip of the STM is contacted to a sample made of 
the same material. The MCBJ experiments start from a 
notched wire that is glued to a bending substrate which 
bends using a piezo-element. At some point the wire 
breaks at the notch forming two fresh surfaces which are 
brought back together to form the atomic contacts. All 
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FIG. 1: Traces of conductance for Ir atomic contacts at 
4.2K showing a jump to contact (left) or a smooth transition 
(right). Arrows indicate the direction of the displacement of 
the electrodes. When there is a jump to contact, the jump 
starts from G a and ends at Gb. 
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experiments were performed at low temperatures (4.2K) 
and under cryogenic vacuum conditions. The traces of 
conductance were taken at constant bias voltage (typi- 
cally at lOmV) which has been shown [9[ to have negligi- 
ble effect on the JC. A current amplifier in the range 
of mA was used to measure the current between the 
electrodes while their distance was controlled using the 
piezo element. Calibration of inter-electrode displace- 
ments were done by means described elsewhere [Tq|. 
Typical traces are shown in Fig. [TJ The conductance 
increases exponentially as the electrodes are brought to- 
gether, followed by a plateau of conductance that re- 
veals the formation of a one- atom contact. In previous 
works 0, i, i, H, 0, 0, H, 0, EH this process has been stud- 
ied experimentally mainly by looking at individual traces 
of conductance. Although this provides some informa- 
tion on the process, it is not suitable for quantitative 
conclusions. Here we have used another approach: First, 
we take a conductance histogram from a few thousands 
of traces of conductance. The conductance histogram 
is used to identify the conductance at which the con- 
tact is formed (G c ). Next, for each contact closing trace 
we automatically record two values, namely the value 
of conductance below G c at which the steepest jump of 
conductance occurs (G a ) and the value of conductance 
immediately following the jump (G5). Finally, these val- 
ues are used to construct a density plot such as the one in 
the left panel in Fig. O The density of pairs (G a ,Gfr) are 
plotted in color coding with G a and G^ along the horizon- 
tal and vertical axis, respectively. The pairs (G a ,G&) for 
which a JC occurs with maximal probability are shown 
as a maximum in this density plot. 

From all metals studied, three characteristic behaviors 
have been identified. Figure [2] shows our results for the 
case of Au from an analysis of more than 300000 traces. 
First of all we notice that there are no data near the line 
G a =Gb, showing that there is always a JC. We observe 
a first maximum at the conductance (0.03,0.87) in units 
of Go, which means that in most of the cases the tun- 
neling regime ends at a conductance of 0.03 Go and then 
the conductance jumps to a value of 0.87 Go- The dis- 
tribution for the jumps is quite broad for this maximum, 
especially in G a . Further analysis shows that this peak 
is most likely the superposition of one at (0.01,0.83) and 
another at (0.05,0.94). We notice a second maximum at 
(0.09,1.6), which corresponds to 7% of the data, as cal- 
culated from the distribution shown at the right panel of 
Fig.0 

In the case of W our observations coincide with those 
of Ref. 6, i.e., the JC is nearly always absent, as seen 
in the left panel of Fig. 03 The evolution of the traces 
of conductance as the contact is stretched is completely 
different from those for Au since there is not a sharp 
transition between conductance plateaus. In the density 
plot we observe most of the data concentrated near the 
line defined by G^ =G a . Although there are some abrupt 
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FIG. 2: (color online) Analysis of the steepest jump of conduc- 
tance before the formation of a metallic contact for the case 
of gold, made from more than 300 000 conductance traces. 
The left panel shows a density plot, where the horizontal axes 
represents the conductance at which the jump takes place 
and the vertical axes shows the conductance of the contact 
formed. We have artificially changed the colors of the peak 
above (gray scale) to make it visible. The right panel shows 
the corresponding histogram of the conductance of the con- 
tact formed after the jump. 



transitions, these show no reproducibility and are seen in 
the density plot as the scarce data away from the G^ =G a 
line. 

Finally, in the case of Ir, shown in the right panel of Fig. [3] 
we find two maxima, one clearly shown at (1.1,2.4) Go 
and another on the G a =G& line. The first maximum 
corresponds to a JC as observed for Au. The conduc- 
tance at which the contact is formed is more than twice 
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FIG. 3: (color online) Analysis of the steepest jump of conduc- 
tance before the formation of a metallic contact for the case 
of W (left), and Ir (right) made from more than 10 000 and 
70 000 conductance traces, respectively. The horizontal axes 
represents the conductance at which the jump takes place and 
the vertical axes shows the conductance of the contact formed. 
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that of Au because of the higher valence of Ir with re- 
spect to Au[l2[. The valence of the metal determines the 
number of electronic channels through the contact, and 
each channel contributes a conductance with a maximum 
of the quantum unit, Go- Note that the most probable 
conductance in the tunneling regime just before the Ir 
contact is formed is higher than one quantum of con- 
ductance. The second peak is located close to the line 
Gb =G a and shows the 21% of cases when the JC is ab- 
sent or is too small to be observed, which is similar to 
the behavior we observed for W. 

These three examples are representative of what happens 
during contact formation. We have analyzed other ma- 
terials and summarize our results in the following table: 



Metal 


G C (G ) 


JC 


maxima (G o ,G 6 )(G ) 


Au 


1.0 


100% 


(0.01,0.87);(0.05,0.94);(0.09,1.6) 


Pt 


1.5 


100% 


(0.13,1.4);(0.19,2.0);(0.35,3.1) 


Ag 


1.0 


100% 


(0.07,1.0);(0.2,1.6) 


Ir 


2.0 


79% 


(1.1,2.4) 


Ni 


1.3 


75% 


(0.2,1.2) 


W 


1.0 


0% 





This table shows, for the different materials the conduc- 
tance of the first plateau, the percentage of cases having 
a JC, and the positions of the maxima that can be clearly 
identified from the density plots. Note that for Pt three 
peaks are observed, like for Au, while for Ag only two 
were resolved. 

The process of JC is related to the mechanical stiffness 
of the contact in competition with the attractive bind- 
ing forces between the atoms to be contacted. Atomistic 
simulations [2, 0, [3] have shown that when surfaces are 
brought into contact a jump occurs due to the adhesion 
forces between the metallic layers. However, this jump 
depends on the stiffness of the material and the geome- 
try of the surfaces involved. For example, we have per- 
formed simple static calculations using embedded atom 
potentials for different surface geometries. The surfaces 
are shifted rigidly each step and then relaxed using con- 
jugate gradient algorithms. From these calculations we 
observe a clear jump between two perfect (100) surfaces 
as well as for an atomically sharp tip approaching a per- 
fect surface. However, the transition is smooth when two 
tips approach each other. 

In the case of the nanocontacts described above, it is 
therefore expected that the existence of a JC will depend 
on the shape of the surfaces before contact. In order to 
understand this effect we have performed molecular dy- 
namics simulations of the elongation of nanowires of Au 
and W until fracture. These fractured structures are then 
dynamically brought into contact. The simulations con- 
sisted of 525 and 437 atoms for Au and W, respectively, 
where three atomic layers at the top and bottom of the 
simulation box were displaced every 1000 time steps by 
a fixed distance, resulting in deformation rates charac- 
teristic of these types of calculations [15]. The simula- 



tion temperature was fixed to 4.2 K. An EAM potential 
was used for the case of Au llq while a Finnis- Sinclair 
potential was used for W [17]. A total of 25 cases were 
computed both for Au and W. The atoms that first come 
into contact were identified considering a cut-off distance 
between first and second nearest neighbors. The distance 
between these atoms is then calculated as the two sides 
approach. Fig. [4] shows the distance between the pairs 
of atoms that first come into contact, for two different 
configurations. Notice that in both cases there is a jump 
at contact. This occurs for all simulations performed for 
Au. The insets in Fig. [4] show the configurations at the 
point of contact for these two cases. In the left panel the 
contact is made by two atoms, one on each side, forming 
a dimer, while in the right panel a single atom (monomer) 
from one side contacts three atoms on the opposite side. 
From the 25 cases simulated 72% dimers are formed at 
contact, 20% monomers, while 8% are contacts through 
two parallel dimers or a double-bond. The existence of 
three different configurations is in agreement with the 
experimental observations for Au, where three peaks are 
observed in the density plot. Moreover, comparing the 
statistics obtained from these calculations with the ex- 
perimental values we infer that the third peak for Au in 
the table above corresponds to a double-bonded configu- 
ration, while the first two peaks correspond to monomer 
and dimer configurations. 
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FIG. 4: Distances between the pairs of atoms that first come 
into contact during the approach for two different configura- 
tions for Au obtained from molecular dynamics simulations. 
The most common configuration (a) consists of a dimer, while 
(b) shows the case of a monomer. The insets show the posi- 
tions of the atoms. Dark atoms are those at contact. 

In order to test the validity of the interpretations 
of the experimental observations, the conductance 
was computed using the first-principles transport code 
ALACANT[18| for representative cases in all three con- 
figurations mentioned above. We have used either local 
density or generalized gradient corrected approximations 
for the exchange-correlation potential and pseudopoten- 
tial basis sets with system sizes of up to 100 atoms. Al- 
though a statistical analysis is not feasible, the results 
lie within the experimental distribution and confirm the 
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above established correspondence between preferred con- 
ductance values and specific atomic configurations, that 
is, the third peak for Au corresponds to a double-bonded 
configuration, while the first two peaks would correspond 
to monomer and dimer configurations respectively. 
Finally, W was studied by molecular dynamics in order 
to understand those cases where no JC occurs. Figure [5] 
shows a particular configuration for W at the point of 
contact. Unlike Au, for W there are configurations for 
which a jump is not observed at contact. These configu- 
rations are closer to two ideal tips than those found for 
Au. For a brittle material such as W much sharper and 
ordered structures are expected than for more ductile ma- 
terials such as Au. In fact, as pointed out previously [6[, 
the elastic properties of these metals play an important 
role in the presence of a JC. 
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FIG. 5: Distances between the two atoms that first come 
into contact for W, as obtained from a molecular dynamics 
simulation. The inset shows the positions of the atoms. Dark 
spheres are those atoms at contact. Contact (defined with the 
cut-off between first and second nearest neighbors) occurs at 
130ps. Notice the different scale with respect to Fig. [4] 

In conclusion, a systematic analysis of the occurrence of 
a jump to contact in different metals has been performed 
using MCBJ and STM experiments. These experiments 
demonstrate the absence of a jump to contact for several 
metals. The results are interpreted using atomistic cal- 
culations of Au and W. The geometry of the electrodes 
before contact seem to play a very significant role deter- 
mining whether a JC occurs, and this shape is related to 
the elastic properties of the metals. Consequently, softer 
metals such as Au or Ag always present a JC while more 
brittle ones such as W do not. Moreover, with the help of 
first-principles transport calculations we have been able 
to correlate the preferential values of conductance mea- 
sured at JC in Au with three distinct configurations: a 
dimer, a monomer and a double bond. 
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